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Abstract 

Given the need for parallel increases in food and energy production from crops in the context of global change, crop 
simulation models and data sets to feed these models with photosynthesis and respiration parameters are increasingly 
important. This study provides information on photosynthesis and respiration for three energy crops (sunflower, kenaf, 
and cynara), reviews relevant information for five other crops (wheat, barley, cotton, tobacco, and grape), and assesses 
how conserved photosynthesis parameters are among crops. Using large data sets and optimization techniques, the C 3 
leaf photosynthesis model of Farquhar, von Caemmerer, and Berry (FvCB) and an empirical night respiration model for 
tested energy crops accounting for effects of temperature and leaf nitrogen were parameterized. Instead of the common 
approach of using information on net photosynthesis response to C0 2 at the stomatal cavity (A n -Cj), the model was 
parameterized by analysing the photosynthesis response to incident light intensity (A n -/ jnc ). Convincing evidence is 
provided that the maximum Rubisco carboxylation rate or the maximum electron transport rate was very similar whether 
derived from A n -Ci or from j4 n -/inc data sets. Parameters characterizing Rubisco limitation, electron transport limitation, 
the degree to which light inhibits leaf respiration, night respiration, and the minimum leaf nitrogen required for 
photosynthesis were then determined. Model predictions were validated against independent sets. Only a few FvCB 
parameters were conserved among crop species, thus species-specific FvCB model parameters are needed for crop 
modelling. Therefore, information from readily available but underexplored 4 n -/ inc data should be re-analysed, thereby 
expanding the potential of combining classical photosynthetic data and the biochemical model. 

Key words: A-I inc curves, acclimation, bioenergy crops, crop modelling, day and night respiration, electron transport rate, leaf 
nitrogen, photosynthesis, Rubisco carboxylation, temperature. 



Introduction 

In conventional crop modelling leaf photosynthesis is 
calculated from net photosynthesis light response curves 
(A n -I inc ; see symbols explanation in Table 1) at ambient 
atmospheric C0 2 level using empirical functions (e.g. 
SUCROS; Goudriaan and van Laar, 1994). In the context 
of better understanding biological processes and exploring 
the impact of climate change, recent crop models (e.g. 
GECROS; Yin and van Laar, 2005), 3D models (e.g. 
Evers et aL, 2010), or terrestrial ecosystem models (e.g. 
LPJmL; Beringer et aL, 2011) calculate photosynth- 
esis based on the mechanistic model of Farquhar, von 



Caemmerer, and Berry (Farquhar et aL, 1980; the FvCB 
model hereafter). 

The FvCB model describes photosynthesis as the mini- 
mum of the Rubisco-limited rate and the electron transport- 
limited rate. The key parameters of the model are the 
maximum Rubisco carboxylation rate (F cmax ), the maxi- 
mum electron transport rate (/ max ), and the mitochondrial 
day respiration (R d ). These biochemical parameters are 
influenced both by the physiological status of a leaf such as 
the amount of leaf nitrogen per unit area (7V a ) (e.g. Harley 
et aL, 1992) and by short- and long-term changes of 
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Table 1. List of main symbols used in this study with their 
definitions and units 



Symbol Definition 



Unit 



A, 
A 

^n,max 

a R 
b R 
Co 
Ci 

E K mc , ^K mo 

^ R n(a) 
E Rn(b, 



J 

^max 

>^max25 

^mc 

^mo 

A/ a 

A/ b 

O 

R 

fin 



e 

K2LL 
^C02LL 

Xj 
Xr 

Xv 



Rubisco-limited net 
photosynthetic rate 
Net assimilation rate 
Electron transport-limited 
net photosynthetic rate 
Light- saturated A n 
x-axis intercept in Equation 12 
Slope parameter in Equation 12 
C0 2 chloroplast partial pressure 
Intercellular C0 2 partial pressure 
Deactivation energy of J max and 



( (Equation 6) 



Activation energy for K mc and for K mo 

Activation energy of J max , R n , and 

1/cmax (Equations 5-6) 

Constant parameter (Equation 1 1) 

Slope parameter in Equation 1 1 

Mesophyll conductance for 

C0 2 diffusion 

Stomatal conductance for H 2 0 

Incident light on leaf surface 

Photosystem II electron transport rate 

Maximum electron transport rate 

Value of J max at 25 °C 

Michaelis-Menten constant for C0 2 

Michaelis-Menten constant for 0 2 

Leaf nitrogen per unit area 

Minimum N a required for 

photosynthesis 

Oxygen partial pressure 

of the air (=210) 

Universal gas constant (=8.314) 

Day respiration rate 

Night respiration rate 

Value of fi n at 25 °C 

Entropy term for J max and V cmax 

(Equation 6) 

Maximum carboxylation rate 
Value of V cmax at 25 °C 
Ci-based C0 2 compensation 
point in the absence of R 6 
Convexity factor for the 
response of J to / inc 
Conversion efficiency of / inc 
into J at low light 
Apparent quantum yield 
of A n at low / inc 
Slope of the J maX 25 and N a 
relationship (Equation 10) 
Slope of the f? n25 and N a 
relationship (Equation 8) 
Slope of the V cmax25 and N a 
relationship (Equation 9) 



nmol C0 2 m 2 s 



nmol C0 2 m s 
limol C0 2 m 2 s 1 



limol C0 2 m 2 s 1 

l^bar 
l^bar 
J mol 1 

J mor 1 
J mol 1 

J mor 1 

J m 2 mol 1 g" 1 N 
mol rrT 2 s~ 1 

mol rrr 2 s~ 1 

lamol photons rrT 2 s _1 

lamol e~ rrr 2 s~ 1 



nbar 
mbar 

g N rrr 2 leaf 
g N m 2 leaf 

mbar 

J K 1 mol 1 
lamol C0 2 rrr 2 s~ 1 
lamol C0 2 rrr 2 s 1 
lamol C0 2 m 2 s 1 
J K~ 1 mor 1 

limol C0 2 rrr 2 s~ 1 
lamol C0 2 rrr 2 s~ 1 
nbar 



mol e~ mor photons 
mol C0 2 mol 1 photons 
limol e" g~ 1 N s~ 1 
\imo\ C0 2 g~ 1 N s~ 1 



l^mol C0 2 g 1 N s 1 



environmental variables such as temperature, light (e.g. 
Hikosaka, 2005), C0 2 (e.g. Makino et al, 1994), and 
drought (e.g. Galmes et al, 2007). 



Usually, the FvCB parameters are obtained by analysis of 
net photosynthesis response to C0 2 at the stomatal cavity 
(A n -Q) (e.g. Sharkey et al, 2007) or by combining A n -Q 
and A n -I inc curves (e.g. Braune et al, 2009) or by combining 
these curves with chlorophyll fluorescence measurements 
(Yin et al, 2009). Obviously, to parameterize the FvCB 
model, information on A n -C[ is predominantly considered to 
be essential, and an ongoing discussion is mainly focused on 
improving the methods of analysing these A n -Q curves 
(Ethier et al, 2004; Sharkey et al, 2007; Gu et al, 2010). 

In the context of forward crop modelling typically for 
predictions at the ambient C0 2 level, the FvCB model is used 
to project leaf photosynthetic rates in response to both 
temporal (diurnal and seasonal) and spatial (within a crop 
canopy) variation in light intensity. This implies that in the 
context of inverse modelling, important FvCB model param- 
eters J max and F cmax should and can be estimated from A n 
responses to I mc . This would reflect better the tradition 
whereby crop modellers describe leaf photosynthesis from its 
response to light intensity (e.g. Goudriaan, 1979), in contrast 
to the tradition that photosynthesis physiologists study gas 
exchange measurements mainly across various levels of C0 2 
(e.g. von Caemmerer and Farquhar, 1981). In fact, the FvCB 
model can be parameterized from analysis of A n -L mc data 
alone (Niinemets and Tenhunen, 1997; Kosugi et al, 2003), 
but so far there is no information about the accuracy of J max 
and F cmax parameters derived from such an analysis. If / max 
and F cmax estimates derived from analysis of A n -I inc are 
similar to those obtained from the common A n -Q analysis or 
combined analysis of A n -Q and A n -L mc curves, it may 
generate an opportunity to reduce empiricism in crop models 
by using readily available A n -I inc data. Therefore, the first 
objective of this study is to explore this opportunity by 
parameterizing the FvCB model using A n -I inc data. 

In the light of current trends for a parallel increase in food 
and energy production from crop species in the context of 
climate change, the use of the FvCB-based simulation models 
together with an urgent need to feed these models with 
photosynthetic and respiration parameters has been increased 
(e.g. Beringer et al, 2011). Compared with the rich in- 
formation found for trees in the literature, there are only 
a few reports on </ max , F cmax , R d , and night respiration (R n ) 
parameters in relation to environmental and management 
factors for economically important crop species (e.g. Muller 
et al, 2005; Braune et al, 2009; Yin et al, 2009) and these 
are virtually lacking for new bioenergy species. Therefore, the 
second objective of this study is 3-fold: (i) to provide new 
information on photosynthesis and respiration for three 
Mediterranean energy crops (Helianthus annuus, sunflower; 
Hibiscus cannabinus, kenaf; and Cynara cardunculus, cynara); 
(ii) to summarize existing information for five major cash 
crops (wheat, barley, cotton, tobacco, and grape); and (iii) to 
assess how conserved FvCB parameters are among crop 
species to better assist modellers in this exploitation. 

Sunflower, kenaf, and cynara crops were chosen because 
these crops have great potential to increase bioenergy 
production in the Mediterranean region (Archontoulis et al , 
2010a, b; Danalatos and Archontoulis, 2010). In addition, 
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the chosen crops cover a wide range of bioindustrial 
applications (biodiesel, bioethanol, heat, and electricity) and 
fit into different cropping strategies (short or long growing 
period, cultivation with or without irrigation, etc.). Sunflower 
is widely grown in the Mediterranean region, but kenaf and 
cynara cultivation is still in the experimental phase; relevant 
information for crop modelling is currently being accumu- 
lated for these crops, including vertical distribution of light 
and nitrogen within crop canopies (Archontoulis et aL, 
2011). Photosynthetic gas exchange studies for sunflower 
have been reported (e.g. Connor et aL, 1993), but there are 
only a few for kenaf (Muchow, 1990; Cosentino et aL, 2004) 
and none for cynara. 

The present analysis focuses on the FvCB parameters in 
response to temperature and 7V a for these bioenergy crops. 
This is because earlier studies on F cmax and J max temperature 
dependencies showed great species-to-species variability 
(Leuning, 2002; Medlyn et aL, 2002(2), and because N a is 
linearly related to Rubisco content that drives C0 2 fixation 
(Makino et aL, 1994), reflects leaf dynamics well (leaf age, 
rank; Archontoulis et aL, 2011), and comprises a reference 
index for scaling photosynthetic C0 2 assimilation from leaf 
to canopy levels (de Pury and Farquhar, 1997). Among 
bioenergy crops, the perennial cynara has long annual 
growth cycles (—10 months each; Archontoulis et aL, 
2010a). Given the numerous reports together with their 
diverse findings on photosynthetic and respiratory acclima- 
tion to growth environment (Atkin et aL, 2005; Ow et aL, 
2008; Yamori et aL, 2005, 2010; Silim et aL, 2010), seasonal 
acclimation effects on photosynthesis and respiration for the 
cynara crop are also investigated. 

Materials and methods 

Literature data for A n -C/ versus A n -\ inc curves 

The first objective of this study was to compare V cmax and J max 
estimates derived either from A n -Q or from A n -I inc curves. For 
this, published data from Yin et al. (2009) for Triticum aestivum 
(cv. Minaret) were used. All relevant parameter values required to 
fit the FvCB model to the A n -Q or A n -I mc data set were available, 
therefore avoiding any statistical artefact in K cmax and / max 
estimation. Wheat measurements (four replicates; all at 25 °C) 
were conducted on leaves with different N a status (15 sets of A n -Q 
and 15 of A n -I mc curves), allowing the comparison of J max and 
J^cmax estimates to be made over a wide range of their values. For 
more information about the measurements, see Yin et al. (2009). 

Energy crop species and study site 

Sunflower (cv. Panter), kenaf (cv. Everglades 41), and cynara (cv. 
Biango avorio) crops were grown in different sections of the same 
field (for details, see Archontoulis et al., 2011) in central Greece 
(39°25'43.4" N, 22°05'09.7" E, 105 m asl) for 3 years (2007-2009). 
The site has a Mediterranean climate with cold/wet winters and 
warm/dry summers (Supplementary Fig. SI available at JXB 
online). The soil was loamy, classified as Aquic Xerofluvent, with 
a shallow groundwater table (1.8-2.8 m below the surface during 
May). In general, crops grown at that site produce much higher 
biomass yields than crops grown on dry soils (e.g. Archontoulis 
et al, 2010/?). During summer, sunflower and kenaf crops were 
frequently irrigated at intervals of 4-6 d according to potential 
evapotranspiration (for site-specific calculations, see Danalatos and 



Archontoulis, 2010) while cynara was irrigated only a few times, 
when necessary during May- June but not during November-April 
(see precipitation in Supplementary Fig. SI). 



Gas exchange measurements and experimental protocol 

Leaf gas exchange (GE) measurements were implemented in situ in 
fully expanded leaves using a portable open gas exchange system 
with a 6.25 cm 2 clamp-on leaf chamber (ADC, LCi/LCpro+, 
Bioscientific Ltd, Hoddesdon, UK). C0 2 /H 2 0 exchanged by the 
leaf was measured using an infrared gas analyser in a differential 
mode. The system allowed for an automated microclimate control in 
the leaf chamber. Before each measurement, attached leaves were 
adapted for 10-45 min to chamber conditions, depending on leaf age, 
time of the day, and season. Daytime GE measurements were taken 
within 1-2 d after irrigation and during morning hours to ensure no 
water stress and to avoid midday depression of photosynthesis. 
Night-time GE measurements were initiated 3CM5 min after sunset 
and lasted for 4—5 h each time. 

To parameterize the model, a common experimental protocol 
was applied per species, including four different sets of GE 
measurements. In all sets, C0 2 concentration was kept at 3 80 ±5 
umol mol -1 . The first set aimed to determine the response of net 
photosynthesis (A n ) to incident light (7i nc ). Accordingly, at fixed 
leaf temperature and measured N a , A n was determined in 11 7 inc 
steps (2000, 1500, 1000, 500, 250, 200, 150, 100, 50, 20, and 0 umol 
photons m -2 s _1 ); in total, 76 curves were constructed. Adaptation 
time to each 7 inc level was ~5 min, except for I[ nc =0 where it was 
>10 min; 3-5 replicated A n measurements were taken at each I inc 
step to ensure stability and precision of measurements. Given that 
the examination of steady-state photosynthesis takes considerable 
time and that GE measurements should be done within a limited 
time frame in order to avoid stress conditions (see above), the 
response of A n to leaf temperature (set II) was determined at three 
Tine levels: 450, 900, and 1800 umol photons m -2 s _1 . N a was also 
determined. At each 7 inc , leaf temperature was increased or 
decreased up to 10 °C from the ambient temperature in steps of 
2-4 °C and replicated A n measurements were recorded every 5 min. 
To establish the relationship between net photosynthesis and N a 
(set III), it was necessary to evaluate leaves with as wide an N a 
range as possible. So, in addition to earlier sets, A n measurements 
were done at saturated 7 inc (1600-1800 umol photons m -2 s _1 ) on 
leaves from different insertion heights in the canopy, from different 
growth stages, and from plots with different N status. Per leaf 
(—180 leaves assessed), 5-10 measurements were taken at leaf 
temperature close to the ambient temperature. 

To obtain direct measurements of the mitochondrial respiration 
occurring in the night (R n ), the response of R n to temperature was 
investigated (set IV). Leaf temperature increased or decreased up to 
10 °C from the ambient temperature in small steps of 1-2 °C, and 
replicated R n measurements were recorded every 4 min. Measure- 
ments were done on leaves with (as much as possible) variable N a . 

To validate the models, GE measurements obtained from the 
same genotypes growing in the same site during summer 2005 and 
2006 (set V) were used. Sunflower and kenaf GE measurements 
were collected using similar techniques and time frames to those 
described for sets I-IV. In cynara, a different protocol was 
followed. The external unit that controls chamber microclimate 
was removed to obtain measurements under real ambient con- 
ditions. Measurements were recorded every 4-8 min, while climatic 
variables were continuously changing following 24 diurnal trends, 
thereby providing a data set to assess whether the FvCB model can 
predict A n under real fluctuating field conditions. 

The wide range of measuring temperature used (15^40 °C) 
unavoidably resulted in variation in vapour pressure difference 
(VPD). An effort was made to reduce that variation by keeping 
humidity high at high temperature. In most cases, VPD was 
maintained below 3 kPa to prevent stomatal closure (Bernacchi 
et aL, 2001). Although VPD was sometimes above 3 kPa at the 
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highest temperatures, the stomatal conductance for H 2 0 vapour 
was not less than 0.30 mol m -2 s _1 (as in Yamori et al, 2005). 

All measured A n data were corrected for the C0 2 respired under 
the gasket surface (total 4 mm width; R. Newman, personal 
communication) following the common approach of Pons and 
Welschen (2002). All GE characteristics were re-calculated according 
to von Caemmerer and Farquhar (1981), for example to provide the 
Q values that are required as input to the FvCB model (see below). 
In addition, the number of replications and observations were in- 
creased to reduce the measurement noise, especially when low C0 2 
exchange rates were measured (e.g. respiration). 

The portion of the leaf used for measurements was cut and its 
area was measured with a Li-Cor area meter. The leaf material was 
then weighed after drying at 70 °C to constant weight and its total 
nitrogen concentration was measured using the Kjeldahl method. 
From these measurements, the leaf nitrogen content 7V a (g N m -2 ) 
was calculated. 

Model and its parameterization 

The FvCB model predicts A n (umol C0 2 m -2 s _1 ) as the minimum 
of two processes (see Fig. 1), the Rubisco carboxylation-limited 
rate (A c ) and the RuBP regeneration- or electron transport-limited 
rate (Aj): 



A n = min(^ c ,^j) 



(1) 



Rubisco-limited photosynthesis is calculated as a function of 
maximum carboxylation capacity ( V C1 



umol C0 2 m 2 s *): 



Q + K mc {\ + 0/K mo ) 



Rd 



(2) 



where Q (ubar) and O (mbar) are the intercellular partial pressures 
of C0 2 and 0 2 , respectively, K mc (ubar) and K mo (mbar) are the 
Michaelis-Menten coefficients of Rubisco for C0 2 and 0 2 , re- 



spectively, and T* (ubar) is the C0 2 compensation point in the 
absence of (day respiration in umol C0 2 m -2 s _1 , which 
comprises mitochondrial C0 2 release occurring in the light other 
than photorespiration; von Caemmerer et al, 2009). 

There are various equations to describe the rate of photosynthesis 
when RuBP regeneration is limiting (Farquhar and von Caemmerer, 
1982; Yin et al, 2004). The most widely used form is given by: 



_ m-r t ) 
Aj - 4C, + sr. d 



(3) 



where J (umol e~ m -2 s _1 ) is the photosystem II electron transport 
rate that is used for C0 2 fixation and photorespiration. / is related 
to the amount of incident photosynthetically active irradiance (I[ nc ; 
umol photons m -2 s _1 ) by: 



K2LL^inc + Jm 



(^2LL^inc + Jmax) ~ ^0J max K2Lhhn 



26 



(4) 



where J max (jimol e~ m~ 2 s _1 ) is the maximum electron transport 
rate at saturating light levels, 9 is a dimensionless convexity 
factor for the response of / to 7 inc , and k 2LL (mol e~ mol -1 
photons) is the conversion efficiency of I lnc into / at limiting light 
levels (Yin and Struik, 2009a; Yin et al, 2009). The formulation 
of Equations 2 and 3 assumes infinitive mesophyll conductance 
(g m ) for C0 2 transfer to chloroplasts, so that Q is used as the 
proxy for the chloroplast C0 2 level (C c ). There is increasing 
evidence that g m might be low enough to allow a significant 
drawdown of C c from C x in most species (Warren, 2004; Flexas 
et al, 2008). However based on the available GE data, it was 
risky to evaluate g m (Pons et al, 2009; von Caemmerer et al, 
2009; Yin and Struik, 2009&), hence the forms of Equations 2 and 
3 had to be used, as in most earlier studies (e.g. Medlyn et al, 
2002a; Kosugi et al, 2003). Omitting g m in the analysis means 



30 T 




-5 line (|jmol photons m" 2 s" 1 ) 

Fig. 1. Main panel: typical net photosynthesis light response curve (/\ n -/inc) at ambient C0 2 concentration. Curve regions for the Rubisco 
carboxylation-limited rate (/\ c -limited, Equation 2; solid line) and the electron transport-limited rate (/\ r limited, Equation 3; dotted line) are 
indicated. Usually, /\ c -limitation occurs above 1500 umol photons m -2 s _1 ; however, it is also possible that the entire A n -I lnc curve is 
described as/\j-limited. Inset panel: representative portion of theA n -l\ nc curve used in calculations of the day respiration (R d ), night 
respiration {R n ), and apparent quantum yield (<D C o2ll)- and <3> C0 2ll were calculated from linear regression analysis to open circles 
while the filled circle represents the value of the For details, see the Materials and methods. 
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that an appropriate consideration is needed in choosing values of 
the Rubisco kinetic constants (see below). 

The temperature responses of respiration and of Rubisco kinetic 
properties (K mc and ^ mo ) are described using an Arrhenius 
function (Equation 5) while the temperature responses of K cmax 
and J max were explored using a peaked Arrhenius function 
(Equation 6); both functions were normalized with respect to their 
values at 25 °C: 



X = X 25 exp 



E X (T- 25) 



298i?(r + 273) 



X = X25exp 



E X (T- 25) 



298i?(r + 273) 



( 29SS X -D X \ 
" eX Py 29SR ) 



t (T+m)S x -D x 
~ cx t>y R(T+273) 



(5) 



(6) 



where T is the leaf temperature (°C); X 25 is the value of each 
parameter at 25 °C (R n25 , K mc25 , K mo25 , and «/ max25 ); Ex is 

the activation energy of each parameter (E Kn , £k mc? ^k mo? E w , and 
Eji in J mol -1 ); D x is the deactivation energy for J max and V cmax 
(Dj and D w in J mol -1 ); S x is the entropy term for J max and V cmax 
(Sj, *S V in J K _1 mol -1 ), and R is the universal gas constant 
(=8.314 J K -1 mol -1 ). Given that Equation 5 is a special case of 
Equation 6, F-tests were performed to determine whether Equation 
6 described temperature responses of K cmax and / max significantly 
better than did Equation 5. When Equation 6 was overparame- 
terized, as often observed in the literature (Dreyer et al, 2001; 
Medlyn et al, 2002a), then S x was fixed at 650 J K" 1 mol -1 
(Harley et al, 1992). 

Rubisco kinetic properties are generally assumed constant 
among C 3 species (von Caemmerer et al, 2009). However, values 
of these constants and their temperature dependency reported in 
the literature vary appreciably, so the choice of Rubisco parame- 
ters is a matter of considerable uncertainty (Dreyer et al, 2001). In 
this work, similar to many other reports (e.g. Medlyn et al, 2002a; 
Miiller et al, 2005), Rubisco parameters reported by Bernacchi 
et al (2001) were selected because these values (i) were estimated 
from in vivo measurements without disturbance of the leaf; and (ii) 
were derived using the Q-based FvCB model and hence are 
compatible with the present analysis assuming an infinite g m (see 
above). The parameter values are: ^ mC 25=404.9 ubar; K mo2 5=218A 
mbar; E^ mc =19 430 J mol -1 ; and E Kmo =36 380 J mol -1 (Table 1). 
Furthermore, using these values, the temperature dependence of 
T* was calculated as (Yin et al, 2004): 



T* = 0.5O 



K mc 



exp -3.3801 



5220 



298i?(r + 273) 



(7) 



where the factor 0.5 is mol C0 2 released when Rubisco catalyses 
the reaction with 1 mol 0 2 in photorespiration. The term in the 
brackets was derived using Bernacchi et al (2001) parameters for 
temperature dependence of maximum carboxylation and oxgenation 
rates of Rubisco. 

The basal capacity of R n25 , V cmax25 , and / max2 5 is linearly related 
to N a (Harley et al, 1992; Hirose et al, 1997; Miiller et al, 2005; 
Braune et al, 2009): 



n25 



XR(^a-^Vb) 

= X,(N dL -N h ) 

= Xj(^a-Af b ) 



(8) 
(9) 
(10) 



where %r (umol C0 2 g 1 N s l ), x v (umol C0 2 g 1 N s and Xj 
(umol e" g" 1 N s -1 ) are the slopes for R n25 , V cmax25 , and / max25 , 
respectively, and 7V b (g N m -2 ) is the minimum value of N a at or 
below which A n is zero. In principle, this N b is practically impossible 
to measure and its estimation depends on the statistical methods 
used and on the available data sets. For instance, different N h 
estimates were found when different data sets were examined (A n or 



K cmax , or J max ; e.g. Harley et al, 1992; Miiller et al, 2005; 
Supplementary Table SI at JXB online) or when 7V b was estimated 
simultaneously with other parameters in optimization procedures or 
when different equations (linear or non-linear) were applied to the 
same data set (Niinemets and Tenhunen, 1997). Given the simplicity 
required in modelling and the lack of biological interpretation of 
different N h values for the same species, a unique N h value (per 
species) was determined beforehand from direct assessments of A n - 
N a plots. Then this estimate was used as input parameter. 

There is some evidence that the activation energy for respiration 
(E^ ) depends on the position of the leaf in the canopy (Bolstad 
et al, 1999; Griffin et al, 2002) and perhaps E Kn is also associated 
with N a since a close relationship between leaf canopy position and 
N a usually exists (Archontoulis et al, 2011). This was tested by 
assuming a linear relationship between E K and N a : 



ER n 



(11) 



and it was checked whether the slope parameter E$> differed 
significantly from zero. 

So far, temperature and nitrogen relationships for R n have been 
described, as extensive GE measurements during the night period 
were available. However, the FvCB requires estimates for R d , 
which is much more difficult to measure. To estimate R d , 
regression analysis was applied to the linear sections of the A n -I inc 
curves for each species (Fig. 1, inset; Kok method; Sharp et al, 
1984). From this analysis, R& was calculated as the j-axis intercept 
of the linear regression and the corresponding R n was estimated as 
the mean of the A n values at 0 umol photons m -2 s _1 . Addition- 
ally, the apparent quantum efficiency at limiting light (<Dco2ll, 
mol C0 2 mol -1 photons) on the incident light basis was calculated 
from the slope of the regression. The I inc range for this regression 
analysis was typically 20-150 umol m -2 s _1 (Fig. 1, inset), while in 
a few cases the I mc range was slightly different, especially for data 
sets obtained at high temperatures. The estimated R d was then 
related to R n as: 



R d = b R x (R n - a R ) 



(12) 



where & R and a R are the slope and the x-axis intercept of the linear 
model, respectively. By assuming that activation energies for R d 
and R n are similar and taking into account the precise quantifica- 
tion of R n based on a large data set, the temperature and nitrogen 
dependencies of R d can be calculated from combining Equations 5, 
8, 11, and 12. This approach allows R d values to be estimated for 
sets II and III (see above) where 7 inc exceeds 350 umol m~ 2 s _1 , for 
which it was not possible to use the Kok method for estimating R d . 



Summary of parameters and statistics 

The basic equations of the FvCB model, Equations 1-4, capture the 
response of A n to Q and to I mc . Coupled with auxiliary temperature 
(Equations 5-7) and nitrogen (Equations 8-12) equations, the model 
also quantifies leaf photosynthesis and respiration (R d and R n ) in 
response to these environmental variables. Data from sets I-IV were 
analysed using step-wise optimization procedures. Per crop, 16 
parameters were estimated following the order: step 1, N h ; step 2, 
Xr> ^R n(a) , ^R n(b) ; step 3, b K , a K ; step 4, 

Xv? Ey, D v , S w ; step 5, k 2LL ; 
step 6, Xj, Ej, Dj, Sj, and 9 (see the Results). Inputs to the model are: 
Q, Tine, leaf temperature, and N a . So, just like using 
A n -Q curves, using A n -I inc data to calculate FvCB model parame- 
ters (e.g. V cmax ) also requires Q as an input to the model, meaning 
that any (short-term) change in stomatal aperture during the A n -I inc 
measurements will have been reflected in the values of Q and thus 
have little effect on the calculation of the FvCB parameters. For 
a similar reason, the direction of changing I mc levels for measuring 
A-l inc curves will also have little impact on parameter estimation 
(see Yin et al, 2011). 

For each step, regression fitting was carried out using the 
GAUSS method in PROC NLIN of SAS (SAS Institute Inc.). To 
investigate seasonal effects of acclimation on photosynthesis and 
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respiration rates of cynara, data sets were split into two periods: 
a cold period with low light from November to April and a warm 
period with high light from May to June (Supplementary Fig. SI 
at JXB online). Then, dummy variables (Zl=l and Z2=0 for 
warm and Z1=0 and Z2=l for cold periods, respectively) were 
introduced into the regression analysis to separate for the effects. 
A dummy variable was also used to best estimate the N h parameter 
(see the Results). 

The goodness of model fit was assessed by calculating r 2 and the 
relative mean root square error (rRMSE). A sensitivity analysis 
was also performed. Model predictions were validated against 
independent data sets (set V). 

Results 

V C max and J max estimates from A n -Cj and/or A n -\ inc 
curves 

Pcmax and J max were estimated for wheat, from either A n -I inc 
or A n -Q curves alone or from the combined data of the two 
curves. The following parameters were set as inputs to the 
model (see Equations 1-4 and 7): K mC 25 and K mo2 s from 
Bernacchi et al (2001); and R^is, ^n25 ? k 2L l ? and 0 per set of 
data from Yin et al (2009). F cmax and J max were successfully 
estimated simultaneously in 40 out of the 45 cases (15 setsx3 
methodologies). In five cases, it was not possible to estimate 
Pcmax from A n -I[ nc curves because in these cases the entire 
curve was limited (Fig. 1). Then we first calculated F cmax 
directly from Equation 2 with observed Q as input for simple 
substitution using data points where 7 inc >1500 umol pho- 
tons m~ 2 s -1 , and secondly by setting F cmax as an input to 
the model, the / max parameter was estimated again. To be 
consistent, results for all A n -I inc curves were presented 
following the two-step approach, because estimates from 
both approaches were very close. 

Figure 2 illustrates F cmax and J max estimates from A n -Q 
and from A n -I inc curves versus the combination of those 
curves. As expected, F cmax and / max estimates obtained from 
A n -Q curves were almost identical to the estimates based on 
the combined data (r 2 =0.97-0.99). However, it was found 
that A n -I inc curves alone also provided sufficient estimates 
(r 2 =0.9 1-0.93) and thus can be considered as an alternative 
to predominant A n -Q curves to parameterize the FvCB 




model. In fact, regression lines in Fig. 2 were matching across 
a very wide range of F cmax and J max values. Even in cases 
where photosynthetic responses to light were entirely Aj 
limited (Fig. 1), F cmax estimates obtained from either A n -I inc 
or A n -Q data were close (Fig. 2). The slight discrepancy of 
the estimates at high F cmax and J max values (Fig. 2) caused 
a lower r 2 for the A n -L mc compared with the A n -Q estimates. 

Step-wise estimation of model parameters for 
bioenergy crops 

Step 1: M b estimation 

Measured light-saturated A n (v4 n max ) responded non-linearly to 
increasing N a in all tested crops (Fig. 3; r 2 > 0.81; P < 0.001). 
An effect of temperature was detected in this relationship only 
at high 7V a (Fig. 3). To estimate the TVb value properly from 
these plots a dummy variables approach was used, in order to 
obtain a unique N h estimate per crop, while allowing the 
equation to vary with different temperatures (optimum versus 
non-optimum temperature ranges; Fig. 3). Derived parameters 
are listed in Table 2. TVb values for all crops were close to 0.4 g 
N m~ 2 , while the lack of N a data below 0.7 g m~ 2 caused 
a high standard error of the N h estimate (Table 2). 

Step 2: R n in relation to temperature and N a 

By combining Equations 5, 8, and 11, R n parameters were 
estimated (Table 3). In cynara, an additional seasonal 
effect was found, with significantly higher R n rates for the 
winter/cold- compared with the summer/warm-growing 
leaves (Fig. 4). Incorporation of this effect into the model 
improved r 2 from 0.68 to 0.72. Of the two R n parameters, 
temperature sensitivity (£ Rn ) was significantly (P < 0.01) 
affected by season, but the slope of the R n -N a relationship 
(Xr) was not (P=0.263); thus, a common Xr value was 
calculated (Table 3). The R n models' goodness of fit was 
satisfactory (r 2 > 0.72; rRMSE < 0.28 across species). 

Step 3: relationship between R d and R n 

Plotting R d versus gave a good linear relationship with 
no significant differences among species (P=0.225; Fig. 5; 
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Vcmax f rom combined A n - C, and A n - l inc curves J max from combined A n - Cj and A n - l inc curves 

Fig. 2. Relationships between V cmax (umol C0 2 m -2 s~ 1 ) and J max (umol e~ m~ 2 s _1 ) estimated from photosynthetic light response 
curves at ambient C0 2 concentration (open circles; A n -I\ nc ) or from photosynthetic C0 2 response curves at saturated light (filled circles; 
/A n -C|) versus estimates obtained from an analysis of combined A n -I inc and A n -Q curves. Data for/\ n -Cj and A n -I inc measurements are 
from Yin et al. (2009) for Triticum aestivum (r/=15). 
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Fig, 3. Relationships between light-saturated net photosynthesis, A n (/ inc >1500 jimol m" 2 s~ 1 ; CO 2 =380 jimol moP 1 ), and leaf nitrogen 
content, N a . Filled symbols refer to data obtained at temperatures near the optimum temperature for photosynthesis per species (sunflower, 
26-34 °C; kenaf, 27-35 °C; cynara, 23-31 °C) and open symbols refer to data obtained at sub- (open squares) or supra- (open triangles) 
optimum temperature ranges. Each point is an average of 4-10 measurements. Lines are fits from a three-parameter non-linear equation: 
A n =A n>maK {2/[1 +exp(-c(/V a -/V b ))]-1 }, (see Sinclair and Horie 1 989), where Ai.max is the asymptote (maximum value) of the dependent variable; 
c is the parameter determining the steepness of the curve; and N b is the intercept of the x-axis denoting a threshold leaf nitrogen value at or 
below A n equals zero. Estimates of parameters are given in Table 2. Cynara's data points were mostly collected during May-June. 



Table 2. Estimates (SE in parentheses) of the non-linear equation 
used to describe data illustrated in Fig. 3 

^n,max is the maximum net assimilation rate (|imol C0 2 m" 2 s" 1 ) at 
saturated light, maximal leaf nitrogen content, ambient C0 2 concen- 
tration, and at optimum (filled symbols) and non-optimum (open 
symbols) temperature ranges; c is a dimensionless factor determin- 
ing the steepness of the non-linear model; and A/ b is the minimum 
leaf nitrogen content (g N m" 2 ) required for photosynthesis. 



Species 


Symbol 3 


^n,max 


c 




Sunflower 


Filled (26-34 °C) 


36.6 (2.48) 


1.19 (0.195) 


0.387 (0.078) 




Open 


26.4 (1.45) 


1.65 (0.313) 




Kenaf 


Filled (27-35 °C) 


35.8 (2.18) 


1 .29 (0.269) 


0.390 (0.126) 




Open 


29.2 (2.12) 


1 .45 (0.334) 




Cynara 


Filled (22-31 °C) 


36.4 (2.41) 


1.08 (0.191) 


0.416(0.097) 




Open 


23.9 (2.16) 


1 .22 (0.725) 





a Symbols in Fig. 3. 

b The confidence limits for N b are: 0.231-0.541 , 0.139-0.640, and 
0.225-0.608 for sunflower, kenaf, and cynara, respectively (P=0.05). 



Table 3). Analysis showed that mitochondrial respiration 
was inhibited by —28% in the light. The observed x-axis 
intercept (a R =0.39) differed significantly from zero 



(P=0.0039), indicating that R n and R d were not entirely 
proportional (Fig. 5). Additionally, no effect of N a 
(r 2 =0.01; P=0.67) but a significant effect of temperature 
(r 2 =0.18; P=0.008) was found on the R d /R n ratio, showing 
that the ratio approached unity at high temperatures. Simil- 
arly, the R n IA QiSnax ratio — ranging from 7% to 11% across 
bioenergy species — was insensitive to changes in N a (P > 
0.05), but increased significantly with increasing temperature 
(r 2 =0.62; P < 0.01; data not shown). 

Step 4: \J C max in relation to temperature and N a 

The relationships of F cmax to temperature and A^ a were 
quantified by fitting Equations 2 and 5-12 to data obtained 
at high light levels (7 inc ^1500 umol m~ 2 s -1 ) to ensure that 
A n is limited only by Rubisco. All required parameters 
(X v? E v , D y , and S y ) were well estimated. Across species, 
there were small differences in x v (<12%; Table 3), and large 
differences in temperature sensitivities >30 °C (Fig. 6a; 
including other crops). Sunflower temperature sensitivity 
was best described by the peaked Arrhenius equation 
(r 2 =0.736; P < 0.001; Table 3), showing an optimum 
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Table 3. Estimates (SE in parentheses) of parameters used to describe temperature and nitrogen sensitivities of photosynthesis and 
respiration rates in three bioenergy crops 

For cynara, when significant differences between warm and cold seasons were found, two estimates are given. For units see Table 1 . 





Parameter 


Sunflower 


Kenaf 




Cynara-warm a 


Cynara-cold a 




Xr 


0.609 (0.006) 


0.954 (0.015) 






0.775 (0.009) 






^R n(a) 


117 912 (1814) 


100 740 (3250) 




— i u yuu \oo i / ) 




146 956 (4281) 






-23 346 (770) 


-15 743 (1455) 








-26 640 (1858) 




n (night) b 


2492 


1403 






3212 






i 2 


0.799 


0.793 






0.724 




RJR n 


b R 






0.843 (0.040) 










a R 






0.390 (0.107) 








^cmax 


Xv 


73.8 (0.94) 


66.7 (0.92) 






65.2 (0.62) 








53 688 (1631) 


61 812 (1402) 






190 831 (33 853) 






205 638 (355) 


0 






158 486 (30 907) 




s v 


650 c 


0 






550 (108.2) 




J 


Xj 


144.2 (3.4) 


122.1 (1.88) 




100 (0.91) 




92.2 (0.88) 




El 


43 295 (5122) 


28 584 d (1131) 






23 111 (971) 






Di 


125 324 (12 653) 


0 d 






204 489 (218) 






s i 


405 (38.47) 


0 d 






650 c 






K 2 LL 


0.255 (0.018) 


0.278 (0.013) 




0.314(0.014) 




0.419(0.011) 




e 


0.607 (0.027) 


0.627 (0.023) 






0.847 (0.011) 






n (day) b 


1366 


2042 






2334 






r 2 


0.928 


0.909 






0.916 




Ratio 


^max/^cmax 


1.95 


1.83 




1.53 




1.41 




^c/^cmax 


0.0057 


0.0103 






0.0085 





3 Warm period=from early May to end of June; cold period=from November to mid-April; see supplementary Fig. S1 at JXB online. 

b Number of data used in the analysis. 

c Fixed value (see the Materials and methods). 

d Alternatively the following parameters: Ej=28 149, Dj=474 614, and Sj=1482 (with a temperature optimum of 41 .7 °C) gave equal 
temperature sensitivities but values were rejected due to a high standard error of the estimate. 
e Normalized to 25 °C. 

temperature for F cmax at 38.7 °C (calculated from Equation 
Al in the Appendix). For kenaf and cynara no optimum 
temperature was observed within the measurement range tested 
(18^1 °C; Fig. 6a). To explore any acclimation of F cmax to 
growth environments in cynara, the model was allowed to 
estimate different parameters for two contrasting seasons. No 
significant effect of the growing season on x v (65.8 versus 64.3; 
^=0.094) or on E v , D Y , and S Y parameters (P=0.247) was 
found, meaning little seasonal F cmax acclimation. 

Step 5: k 2 ll in relation to temperature and N a 

K21X was estimated indirectly from ®co2ll information (see 
Equation A2). Correlations of k 2LL with temperature, light, 
and nitrogen were investigated afterwards. The results 
indicated poor correlations with N a (r 2 =0.26, P=0.025), leaf 
temperature (r 2 =0.19, P=0.104), and the combination of 
the above (r 2 =0.44, P < 0.01; data not shown). However, 
better relationships were obtained when k 2 ll was regressed 
against seasonal temperature (r 2 =0.40, P=0.004) and radi- 
ation data (r 2 =0.34, P=0.003), showing a long-term k 2 ll 
acclimation. This became clearer when average k 2 ll values 
per crop and per growth environments were considered 
(Fig. 7). These findings were supported fairly well by 
literature data (Fig. 7). Based on this analysis, average k 2LL 
values per species were considered in further analyses 
(including acclimation effect for cynara, Table 3). 



Step 6: J max in relation to temperature and N a 

All J max temperature sensitivities (except kenaf; Table 3) were 
best described using Equation 6. Across species, / max 
temperature sensitivity was highly variable (Fig. 6b including 
other crops), while the maximum / max was obtained at lower 
temperature than the maximum F cmax (temperature optimum 
of 32, 42, and 33 °C for sunflower, kenaf, and cynara, 
respectively; Fig. 6). As a result, there was a decreasing trend 
of the / max /^cmax ratio with increasing temperature (Fig. 8). 
For cynara, a significant (P < 0.05) temporal change was 
found for the Xj parameter (Table 3). Xj showed a larger 
variability (36% change) than x v (12% change) among species 
and growth environments studied (Table 3). The parameter 0 
was lower for sunflower (0.60) and higher for cynara (0.84), 
but close to the commonly used value of 0.75 in all cases. All 
these differences (including temperature and nitrogen sensi- 
tivities) among species and growth environments became 
smaller when the / max /^cmax ratio was plotted against leaf 
temperature (Fig. 8). 

Sensitivity and validation analysis 

To investigate the uncertainty introduced into the estimates 
by the chosen Rubisco kinetic parameters, the initial values 
of Bernacchi et al. (2001) were increased or decreased by 
20% and optimization procedures were repeated. Not 



Temperature and nitrogen effects on leaf photosynthesis and respiration of biomass crops | 903 



A measured during warm periods 
A measured during cold periods 
A predicted 



(a) 




25 30 35 

Leaf temperature (°C) 

Fig, 4. Cynara's night respiration rates (R n ) in relation to leaf 
temperature. Data are presented per growth season and include 
leaves with various A/ a . (a) The predicted R n from a simple tempera- 
ture-sensitive model (Equation 5; parameter values used are shown), 
(b) The predicted R n from a combined nitrogen-, temperature-, and 
acclimation-sensitive model (see parameter values in Table 3). 
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Fig. 5, Relationship between day (f? d ) and night {R n ) respiration 
rates (see also Table 3 and Equation 12). 

surprisingly, a maximum change was obtained in the 
estimated F cmaX 25, whereas the remaining parameters were 
less affected (<5%; data not shown). Given that even the 
maximum change in F cmax was —11% in response to a 20% 
change, the parameter estimates were quite stable despite 
the uncertainties in values of Rubisco kinetic constants. 
A further analysis showed that the predicted A n was sensitive 



to a 20% decrease in % v and Xj, whereas its sensitivity to other 
changes was weak (Fig. 9). 

Lastly the models were validated against independent data 
sets (Fig. 10). Predictions versus observations for sunflower 
and kenaf were satisfactory (rRMSE < 0.15; Fig. 10a, b). For 
cynara the FvCB model was tested using measurements 
from a series of 24 h diurnal cycles (Fig. 10c), where stress 
conditions were unavoidably present (data sets outside the 
calibration range). In general, predictions were close to actual 
measurements, except for those data obtained from 14:00 h 
to 18:00 h, where a systematic overestimation was detected 
(Fig. 10c). The FvCB model responded to lowering temper- 
ature in late afternoon by increasing A n ; however, actual 
measurements indicated that the photosynthetic apparatus 
could not recover so quickly from the 'photosynthesis 
midday depression'. The failure in predicting the depression 
and its after-effect during the recovery hours (Fig. 10c) might 
be attributed to the 'steady-state' character of the FvCB 
model. These results suggest that prediction of diurnal 
photosynthesis for species grown in the Mediterranean region 
requires more detailed approaches in which g m , recovery 
functions for A n (midday depression), and the effects of leaf 
water potential should be included (see Tuzet et al, 2003; 
Vico and Porporato, 2008; Yin and Struik, 2009a). 



Discussion 

Use of A n -\ inc curves to parameterize the FvCB model 

The FvCB model parameters, «/ max and F cmax in particular, 
have been predominantly estimated from A n -Q data sets 
(Harley et al, 1992; Medlyn et al, 2002a). The value of A n -Q 
curves for parameterizing the FvCB model is confirmed (Fig. 2). 
It was also shown that F cmax and J max can be estimated 
sufficiently well by an appropriate analysis of A n -I inc data 
alone (r 2 =0. 9 1-0.93; Fig. 2). Unlike J max , F cmax cannot 
always be estimated from A n -I inc curves; that is when the 
entire curve is Aj limited (Fig. 1). This is often observed in 
field crops (e.g. cotton; Wise et al, 2004). Actually, Boote 
and Pickering (1994) used only the Aj equation of the FvCB 
model to calculate leaf photosynthesis in their canopy 
photosynthesis model. For the purpose of using the 
complete FvCB model, the two-step approach is proposed 
to estimate both F cmax and J max from A n -I inc data. This is 
in line with the approach of Niinemets and Tenhunen 
(1997), but in contrast to that of Kosugi et al (2003) and 
Miiller et al (2005) who assumed a fixed / ma x^cmax ratio of 
2.1 at 25 °C (based on Wullschleger, 1993) in their analyses. 
This assumption does not allow for the flexibility of the 
ratio as observed for different species or for the same 
species when grown under different environments, thereby 
introducing many uncertainties in parameter values (see 
Fig. 8 and discussion below). 

The present results indicated that information from A n -I inc 
curves has been underexplored. Use of A n -I inc curves has an 
additional advantage in that data of A n -Q curves may be 
uncertain due to C0 2 leakage during gas exchange 
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Leaf temperature (°C) Leaf temperature (°C) 

Fig. 6. Temperature sensitivities for V cmax (a) and J max (b). Values normalized to 1 at 25 °C. Filled symbols refer to bioenergy crops while 
open symbols refer to wheat (de Pury and Farquhar, 1997), barley (Braune et al., 2009), cotton (Harley et al., 1992), grapevine (Schultz, 
2003), tobacco (Bernacchi et al., 2001; 2003), and perennial Plantago asiatica (Ishikawa et al., 2007). For kenaf, both observed J max 
temperature sensitivities are plotted (see Table 3; note that for kenaf the measurement range was up to 41 °C). 
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Fig. 7. Conversion efficiency of incident light into linear electron flux (k 2 ll) in relation to (a) seasonal growth temperature, (b) short-term 
changes in leaf temperature, (c) seasonal irradiance, and (d) N a , for three bioenergy crops (filled symbols: see key in Fig. 5) and four 
major field crops (open symbols; see panels for details). Growth temperatures and irradiances were calculated from Supplementary Fig. 
S1 at JXB online. For sunflower and kenaf one average k 2L l value (± vertical standard error) was calculated because measurements 
(set I) were conducted during the July-August period when temperature and radiation do no change much (Supplementary Fig. S1 at 
JXB online). In contrast, for cynara (cold and warm) four average k 2 ll values were calculated, reflecting the months November, April, 
May, and June, respectively. Horizontal bars (when larger than symbols) indicate the mean standard error of the explanatory variable. 
Information on growth temperature and irradiance could not be retrieved from from the studies Muller et al. (2005) and Schultz (2003). 



measurements when C0 2 set point values are either below or 
above the ambient air C0 2 level (Flexas et al, 2007). Crop 
modellers used to measure A n -I xnc curves under an ambient 



C0 2 condition, upon which an empirical model for light- 
response curves is parameterized. Provided that values of Q 
across I inc levels are properly monitored, re-analysing readily 
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Fig, 8. J m ax/^cmax ratio versus leaf temperature. Closed symbols 
refer to bioenergy crops, open symbols to major field crops. Note 
that for cynara two lines were plotted because the parameter 
^max25 differs between seasons (see Equation 10 and Table 3). 



available A n -I inc data to parameterize the FvCB model will 
strengthen photosynthesis calculations in crop models. This 
would expand the potential of combining classical photosyn- 
thetic data and the biochemical FvCB model to assess the 
impact of climate change on crop production and to examine 
options of bioenergy production under a changing climate. 

On the other hand, caution should be exercised that use of 
^n-^inc data sets does not allow the model to account for the 
TPU (triose phosphate utilization)-limited rate, the third 
limitation added by Sharkey (1985) to the FvCB model. TPU 
limitation sets an upper limit to the maximum photosynthetic 
capacity and is usually observed at high CO2 or/and low O2 
levels (e.g. Wise et al, 2004), although many studies still 
ignore this limitation (e.g. Wohlfahrt et al, 1999). The 
limitation, if it occurs, can be easily identified, at the high 
end of A n -Q curves, versus the Rubisco limitation that can 
be identified at the low end of A n -Q curves. In the present 
study, where essentially A n -I inc curves were used, it was not 
possible to detect this limitation, because both Rubisco and 
TPU limitations, if any, will occur at the high end of A n -I inc 
curves. This is certainly the disadvantage of using A n -I inc 
curves to parameterize the FvCB model. Fortunately, the 
present light response curves were obtained under ambient 
C0 2 conditions, so any TPU limitation, if it exists, can be 
assumed to be negligible under these measurement condi- 
tions. In the future climate where the ambient C0 2 level is 
expected to increase, the TPU limitation will be more likely 
to occur. Therefore, use of A n -I mc curves to estimate FvCB 
model parameters needs to be tested across high C0 2 levels 
and a broad range of other environmental variables in order 
to decide how conserved these parameters are. 

Below the effects of temperature, N a , and season on 
photosynthesis and respiration parameters, all derived from 
the current A n -I inc data for three bioenergy crops, are 
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Fig. 9. Sensitivity analysis of the predicted A n in response to 
a ±20% change in input parameter values for the photosynthesis 
model. The relative change in predicted value was calculated as: 
100x(/\ n , predicted-/\ n , predicted, original)//\ n , predicted, original. 
When input parameters were part of a linear or polynomial 
equation (e.g. Ej, Dj, Sf, Equation 6) and strongly intercorrelated, 
a combined change was implemented. 

discussed. The present findings will be compared with those 
reported for the crops wheat, barley, cotton, tobacco, and 
grapevine based on A n -Q or combined A n -Q and A n -I[ nc 
data sets, with attention to any conserved nature in these 
parameters among species. 



Night and day respiration parameters: % R , E R . E Rr 
b R , and a R 



n(bf 



This study is among few in the literature providing direct R n 
measurements, underlining the great importance of respiration 
in carbon budgets (Valentini et al, 2000). The present estimates 
for Xr (range: 0.61-0.95 umol C0 2 g -1 N s -1 ; Table 3) agree 
well with previous reports for crops (Hirose et al, 1997; Reich 
et al, 1998; Miiller et al, 2005; Braune et al, 2009), but 
current values are almost double compared with those for trees 
(Bolstad et al, 1999; Griffin et al, 2002). The temperature 
sensitivity for respiration (E Kn ) was significantly correlated 
with 7V a in all species (Equation 11; Table 3), indicating that 
respiration in leaves with high 7V a values (young/sun leaves) 
was less sensitive to changes in temperature, while leaves with 
lower 7V a values were more sensitive (senescence/shade leaves). 
Griffin et al (2002) and Bolstad et al (1999) working with tree 
leaves that were positioned in different canopy layers — also 
having different 7V a values — found temperature sensitivities 
similar to those in the present study, while Turnbull et al 
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Fig. 10. Measured versus predicted photosynthesis (all panels) and measured versus predicted night respiration (only in the lower 
panel c). In c, canopy C0 2 varied from 350-380 jimol mol -1 during day time to 450-600 jimol mol -1 during night-time; VPD followed 
temperature variations, and stomatal conductance ranged from 0.05 mol rrT 2 s _1 during the night up to 0.48 mol m~ 2 s~ 1 during the 
day time. The model predicted diurnal trends moderately (^=0.814, rRMSE=0.553, n=720). When midday measurements were 
excluded (14:00-16:00 h), the model fit was improved (^=0.930, rRMSE=0.335, n=543). 



(2003) reported the opposite. However, in none of these studies 
was £ Rn significantly correlated with N a . 

For cotton, Harley et al (1992) reported a simple 
temperature-sensitive R n model for leaves with variable 7V a . 
The present analysis indicated that it is useful to calculate 
both R n components as a function of N a (e.g. Fig. 4; across 
all species, r 2 scaled from 0.53 to 0.77). The component R n 25 
accounted for 27% and E Kn for the other 5% of this 
improvement in r 2 . However, the remaining unexplained 
variability in night data sets (see r 2 in Table 3; Fig. 4) means 
that apart from 7V a , other factors should be explored. 

Unlike for R n , it is difficult to measure directly as 
such measurements require sophisticated methodologies (e.g. 
Haupt-Herting et al, 2001; Pinelli and Loreto, 2003; Parnik 
and Keerberg, 2007). Its value is empirically estimated indir- 
ectly using various methods (for a comparison see Yin et al, 
2011), or is commonly fixed as 1% of F cmax or as 50% of R n 
(de Purry and Farquhar, 1997; Wohlfahrt et al, 1998; Medlyn 
et al, 2002a; Kosugi et al, 2003; Braune et al, 2009). Here, 
application of the Kok method (Sharp et al, 1984) indicated 
a 28% reduction in R d compared with R n , an estimate which is 
positioned at the lowest reported range (light inhibition range: 
24-90%; Buckley and Adams, 2011, and references therein). 

Rubisco and electron transport parameters: M b , % v , %j, 
k 2 ll, 9, E w E y , D v , and D y 

The present findings for N h (Fig. 3) along with published 
data support the idea that this threshold value for photosyn- 



thesis is not affected by temperature (Sage and Pearcy, 1987; 
Makino et al, 1994; Niinemets and Tenhunen, 1997), C0 2 
(Harley et al, 1992; Hirose et al, 1997), or irradiance levels 
(Makino et al, 1997). Excluding the statistical bias that 
usually exists in N\y estimations (see the Materials and 
methods) it is believed that a common TVb is 0.3-0.4 g N m~ 2 
for C3 crop species (excluding legume crops; Supplementary 
Table SI at JXB online). For use in modelling, it was shown 
that a ±20% change in the 7V b value resulted in a <5% 
change in the predicted A n (Fig. 9). 

The relationships between A n ^ max and 7V a at near-optimum 
temperature ranges for sunflower, kenaf, and cynara (Fig. 3) 
agreed well with several non-legume C 3 species (Supplementary 
Fig. S2 at JXB online). The observed decline in A n ^ ax at high 
temperature (Fig. 3; 7V a >2 g N m~ 2 ) is associated with g m 
(Bernacchi et al, 2002) and/or F cmax and / max limitations of 
photosynthesis (Fig. 6). 7V a and leaf temperature explained 
>81% of the temporal (seasonal) and spatial (within a crop 
canopy) variation in ^4 n?max values (Fig. 3). The remaining 
unexplained variability might be due to leaf adaptation to 
different microenvironments created by C0 2 and light gra- 
dients within crop stands (Buchmann and Ehlinger, 1998; 
Archontoulis et al, 2011). This may have an additional impact 
on -/max and F cmax estimates and their ratio. 

Nevertheless, the observed consistency among ^4 n ,max-^ a 
plots (Supplementary Fig. S2 at JXB online) along with the 
similar Xv estimates for sunflower, kenaf, cynara, cotton, 
wheat, and barley (range: 60-82 jimol C0 2 g -1 N s -1 ; Table 3; 
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Harley et aL, 1992; de Pury and Farquhar, 1997; Miiller 
et aL, 2005, 2008; Braune et aL, 2009) suggests that x v is 
very conserved for this plant group (v4 n max =30-35 umol 
C0 2 m~ 2 s -1 ; Supplementary Fig. S2). 

Unlike Xv? Xj f° r the same group was highly variable (90- 
165 jimol e~ g -1 N s -1 ). However, the parameter Xj (which 
determines «/ max2 5; Equation 10) is not independent of, but 
interrelated to, the values of k 2L l and 0 (see Equation 4). 
This means that use of constant k 2LL and 9 values across 
species and environments will bias / max estimates and 
therefore the / max /^cmax ratio. Among sunflower, kenaf, 
and cynara, Xj varied by 36%, k 2 ll by 39%, and 0 by 28%, 
but in different directions (Table 3). When k 2LL was fixed to 
0.3 and 0 to 0.7 (commonly assumed values; de Pury and 
Farquhar, 1997; Medlyn et aL, 2002a), the Xj variation 
among crops and growing environments became smaller 
(15%), and the / max /Fcmax ratio less variable. 

The present analysis showed that variation in the electron 
transport rate among bioenergy crops followed changes in 
environmental conditions during growth (Supplementary 
Fig. SI at JXB online), with higher J rates for cynara in low 
light (<700 jimol m~ 2 s -1 ; winter period) and higher J rates 
for sunflower and kenaf in high light conditions (>700 jimol 
m~ 2 s -1 ; summer period; Table 3, Equation 4). This is 
consistent with recent findings for tobacco (Yamori et aL, 
2010) where plants grown under low light enhanced the 
efficiency of light acquisition while those grown under high 
light enhanced the capacity of light utilization, through 
changes in chlorophyll contents, the chlorophyll alb ratio, 
and cytochrome / and Rubisco contents. 

In studies of Wullschleger (1993), Dreyer et aL (2001), 
and Medlyn et aL (2002a) the k 2LL was fixed as a constant 
at 0.18, 0.24, and 0.30, respectively, across species, crop 
stages, and environments. However, Yin et aL (2009) 
directly demonstrated a positive relationship between k 2LL 
and 7V a , which was confirmed by the results of a model 
curve-fitting procedure (Miiller et aL, 2005; Braune et aL, 
2009; Yamori et aL, 2010). 

In Fig. 7, k 2LL information for eight crops is summarized 
and this large variation is interpreted in the light of long- or 
short-term response to temperature or irradiance. Across 
species, the highest k 2LL values were found in crops grown 
under long-term low irradiance and temperature conditions 
(Fig. 7a, c). To understand this, it is necessary to underline 
the components of the k 2LL parameter (see Equation A3 
derived by Yin et aL, 2004, 2009; Yin and Struik, 2009a; 
also see equation 6 in Niinemets and Tenhunen, 1997). The 
fraction of 7 inc absorbed by the leaf photosynthetic pig- 
ments (parameter in Equation A3) is affected by long- 
term changes in light and temperature through its changes 
in leaf morphology. Leaves grown at high temperature are 
generally thinner, with a lower ability to absorb light 
(Poorter and Evans, 1998; Yamori et aL, 2005), therefore 
providing a reasonable explanation for the observed k 2LL 
reduction with increasing temperature. On the other hand, 
leaves grown at high irradiance are thicker (Niinemets and 
Tenhunen, 1997), indicating that k 2LL variation is much 
more complex and still not fully understood. Nonetheless, 



caution should be exercised when modelling canopy photo- 
synthesis based on the sun/shade approach (de Pury and 
Farquhar, 1997; Yin and van Laar, 2005) because k 2LL 
increases with increasing N a (Fig. 7d), while k 2LL also 
increases with decreasing light (Schultz, 2003; shade leaves 
which generally have low N a values; Fig. 7b). 

The normalized temperature functions of F cmax and «/ max 
were variable across crops (Fig. 6), particularly above 
30 °C, in line with Leuning (2002). This mean that the 
assumption used in crop modelling, a unique A n response to 
temperature across crop species, is inappropriate when 
photosynthesis is calculated by the FvCB model. In the case 
of no available data, it is suggested that researchers as a first 
approximation use F cmax and J max temperature parameters 
from species that belong to the same family (see Fig. 6; 
cotton and kenaf belong to Malvaceae; sunflower and 
cynara to Asteraceae). 

The /max'^cmax ratio provides an estimate of the relative 
activities of RuBP regeneration and Rubisco carboxylation, 
and incorporates both temperature and 7V a effects. This 
study confirms (Table 3) the generally reported / max /F cmax 
value of 2.0 ±0.5 (Wullschleger, 1993; Poorter and Evans, 
1998; Bunce, 2000; Leuning, 2002; Medlyn et aL, 2002a). 
However, this ratio should not be considered constant in 
absolute terms. F cmax is dependent on the Rubisco para- 
meters used (up to 11% change; see also Medlyn et aL, 
2002a) and J max is affected by the assumed k 2LL and 0 values 
used (see earlier discussion). For instance, grape showed 
a much higher / max /^cmax ratio compared with other crops 
(Fig. 8). Apart from the effect of species, there are two 
possible artefacts causing this: the different Rubisco para- 
meters used in that study (Schultz, 2003) and the lower grape 
k 2LL values compared with the other crops (Fig. 7b). Also 
use of Q instead of C c affects this ratio. Thus approaches 
(e.g. Kosugi et aL, 2003; Miiller et aL, 2005) that fix the / max / 
^cmax ratio at a constant value to parameterize the FvCB 
model should receive critical reservation. 

Seasonal effects on photosynthesis and respiration in 
cynara 

Direct interpretation of the seasonal effects on A n and R n for 
cynara is difficult because both the climate (Supplementary 
Fig. SI at JXB online) and the plant stage are different, with 
new and old leaves being present (Archontoulis et aL, 2010a; 
Searle et aL, 2011). R n acclimated to cold and warm 
environments to a larger extent than did A n (Table 3; Fig. 4). 
This is consistent with previous studies (Yamori et aL , 2005; 
Ow et aL, 2008; Silim et aL, 2010). 

The nature of R n acclimation is variable within and 
among plant species, and it is usually related to changes in 
E Kn and/or to changes in R n2 5 (Atkin et aL, 2005; Searle 
et aL, 2011). Given that Xr did not change between seasons 
(P=0.269; Table 3) and that the measured winter leaves had 
higher 7V a values than the summer leaves (on average 2.48 
versus 1.53 g N m~ 2 ; see also Fig. 7d), this indicates that 
basal capacity, R n i5, plays an important role in this 
acclimation. Secondly, £ Rn was also higher during winter 
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periods. Apparently, cynara follows an 'acclimation type IF 
(Atkin et al, 2005) where the overall elevation of the R n - 
temperature response was affected by season and growth 
stage (Fig. 4). 

Among FvCB parameters analysed, seasonal effects were 
found on two electron transport parameters, Xj an d k 2 ll 
(Table 3 and earlier discussion), and none related to F cmax . 
Literature information on A n acclimation is diverse among 
studies (Wilson et al , 2000; Medlyn et al , 2002a; Bernacchi 
et al, 2003; Hikosaka, 2005; Yamori et al, 2005; Braune 
et al, 2009; Silim et al, 2010). As far as is known, only 
Wilson et al (2000) reported both Xj an d Xv seasonal 
changes in trees, while Braune et al (2009) found only Xj 
variation for barley as in the present study. For cynara, the 
normalized F cmax and J max temperature functions were 
slightly changed between seasons, in line with other field 
studies (Medlyn et al, 20026; Schultz, 2003), but in contrast 
to growth chamber studies (Bernacchi et al, 2003; Yamori 
et al, 2005; Ishikawa et al, 2007; Braune et al, 2009) where 
plants were grown only at different temperatures. The fact 
that this study assessed leaves with different N a status may 
be a reason, but an inconsistency between actual field and 
controlled chamber studies is obvious. 

The / m ax/^cmax ratio has been reported to be either 
sensitive or insensitive to growth temperature (see discus- 
sion by Hikosaka et al, 2005), growth irradiance (Poorter 
and Evans, 1998; Yamori et al, 2010), and seasonal changes 
(Bunce, 2000; Medlyn et al, 20026). The present results 
suggest that cynara regulates the balance between RuBP 
regeneration and Rubisco carboxylation to maintain the 
4ax/^cmax ratio almost constant (change <8%; Table 3) 
across seasons and growth stages. 



Conclusions 

This study provides new information on photosynthesis and 
respiration rates for three bioenergy crops, sunflower, 
kenaf, and cynara. It provides an alternative way to 
parameterize the FvCB model from A n -I inc data, instead of 
using A n -Q data that are more expensive to obtain. It was 
shown that major FvCB model parameters, F cmax and / max , 
derived from either A n -C[ or A n -I inc analysis, are very close 
(r 2 =0.92). Present models can predict photosynthesis under 
varying levels of Q, 7 inc , temperature, and leaf nitrogen, and 
can estimate night respiration under varying levels of 
temperature and leaf nitrogen, for the three bioenergy 
crops. Comparisons of FvCB model parameters among 
sunflower, kenaf, cynara, cotton, wheat, barley, tobacco, 
and grapevine indicated that only a few parameters were 
conserved. This means that in order to feed crop models 
properly, species-specific FvCB model parameters are 
needed. In this context, readily available A n -I inc data — that 
have been underexplored — can assist in that respect. By 
combining classical photosynthetic data and the biochemi- 
cal model, the potential of crop growth models to assess the 
impact of climate change on crop production and to 
examine options of bioenergy production under a changing 



climate is enlarged. Further research is needed to quantify 
reliably the effects of photosynthetic acclimation and di- 
urnal midday depression identified in this study. 



Supplementary data 

Supplementary data are available at JXB online. 

Figure SI. Average monthly temperatures, radiation, and 
precipitation at the experimental site (period: 2007-2009). 
Sunflower measurements were taken from July to August; 
kenaf measurements from July to September, and cynara 
measurements from November to June. 

Figure S2. Reported relationships between light-saturated 
net assimilation rate at ambient C0 2 concentration and at 
near-optimum temperature (v4 n max in jimol C0 2 m~ 2 s -1 ) 
and leaf nitrogen content (g N m~ 2 ) for C 3 crops (a), C 3 
legume crops and trees (b), and C 4 crops (c). (d) An average 
relationship for C 3 and C 4 crops. 

Table SI. Reported values (minimum leaf nitrogen for 
photosynthesis, in g N m~ 2 ) for various species. 
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Appendix 

1. Estimating the optimum temperature from the peaked 
Arrhenius equation 

The optimum temperature for K cmax or / max in Equation 6 is given 
by the following equation (Medlyn et al, 2002a): 



D x 



1 opt 



S x -R-\n(^ 



(Al) 



2. The relationship between k 2 ll and ®co2/_/_ 
By dividing both parts of Equation 3 by I mc and re-arranging, the 
efficiency of incident light conversion into e~, k 2 ll ? can be 
calculated mathematically from <3> C0 2ll' 



Rn 



C0 2 LL- 



K 2LL~ A 



K2LL = @C0 2 LL 



4Q + 8r» 
4C,- + 8 A 



(A2) 



This approach was also used by Niinemets et al (2001), but 
lacks any further interpretation. 
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3. Components of parameter k 2 ll 

Yin et al. (2004) described a generalized stoichiometric equation 
for Aj, where the linear photosystem II (PSII) electron transport 
rate (J) was replaced by the total electron transport rate passing 
PSII (J 2 ) and fractions of the total e~ flux passing PSI that follow 
cyclic (/ C yc) and pseudocyclic (/p Se udo) pathways. Again, under low 
light conditions, dividing J by I inc yields k 2 ll as follows (Yin and 
Struik, 2009a; Yin et al, 2009): 

-2LL = {- |^=f f 1 - £f~) = P 2 P*2 LL ( 1 - ^f] 
1 inc 1 inc \ 1 J eye I \ 1 J eye I 

(A3) 

By definition, the variable J 2 can be replaced by the term 
P2X/?x<p2LL x ^inc 9 where p 2 is the fraction of absorbed 
irradiance partitioned to PSII (usually assumed to be 0.5), 
is the fraction of I[ nc absorbed by the leaf photosynthetic 
pigments, and 0 2 ll is the PSII e~ transport efficiency under 
limiting light. 
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